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Radiometr ically  and  photogr ammetr ically  calibrated  photo¬ 
graphic  images  of  tracer  material  released  into  the  upper  atmo¬ 
sphere  were  microdensitometrically  reduced.  The  digital  optical 
density  data  were  computer  analyzed  by  Fourier-transform  and 
configuration  space  techniques  to  determine  turbulent  diffusion 
coefficients  and  heating  rates,  and  other  parameters  describing 
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20.  Abstract  (continued) 

upper  atmosphere  turbulence.  Turbulent  release  altitudes 
were  determined  by  photogrammetric  triangulation  and  by  an 
inverse  triangulation  method  using  previously  determined 
wind  profiles.  Examination  of  the  long  term  variation  of 
turbopause  height  suggests  that  the  altitude  of  the  turbo¬ 
pause  is  a  function  of  the  10.7  cm  solar  flux. 
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FOREWORD 


The  program  described  here  consisted  of  quantitative 
characterization  of  turbulence  present  in  90  to  120  kn  altitude 
region  of  the  atmosphere.  This  has  been  accomplished  through 
computer  analysis  of  photometrically  and  photogrammetr ically 
calibrated  images  of  chemical  tracer  releases  to  determine 
Fourier  turbulence  spectra,  configuration  space  statistical 
descriptors  of  turbulence  and  atmospheric  winds.  Data 
acquired  by  in  situ  measurement  techniques  (temperature  and 
wind)  were  often  used  to  supplement  information  derived  photo¬ 
graphically.  Similar  earlier  work  performed  by  the  PhotoMetrics 
research  group  is  described  in  References  1-4,  and  References 
5-10  are  examples  of  the  application  of  the  chemical  and  smoke 
trail  methods  to  measurement  of  atmospheric  turbulence  and 
dynamics . 

The  author  wishes  to  express  his  thanks  to  Mr.  Samuel 
P.  Zimmerman  (Technical  Monitor)  of  AFGL  for  his  continued 
encouragement  and  support  and  to  K.  McDevitt  and  C.  Rice 
for  their  contributions  to  this  report. 
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SECTION  I 


INTRODUCTION  AND  SUMMARY 

The  work  reported  here  was  directed  toward  the  extraction 
of  turbulent  diffusion  coefficients  and  heating  rates  from 
photometrically  calibrated  photographs  of  chemical  smoke  tracers 
released  into  the  atmosphere.  Previously  reduced  data  related  to 
turbulent  diffusion  and  transport  were  also  evaluated;  specifi¬ 
cally  atmospheric  wind  profiles  determined  by  photogrammetr ic 
triangulation  and  temperature  and  wind  profiles  measured  by 
falling  spheres.  Further  effort  was  devoted  to  an  analysis 
of  the  long  term  variation  of  turbopause  height  as  a  function 
of  the  solar  10.7  cm  flux. 

The  surface  brightness  and  transport  of  the  chemical 
tracer  clouds  are  derived  by  microdens i tometr ically  scanning 
photographically  recorded  images  of  their  self-radiation  and 
scattering  of  sunlight.  The  reduced  data  are  analyzed  to 
determine  both  Fourier  transform  and  configuration  space 
descriptors  of  the  fluctuating  component  of  the  surface  bright¬ 
ness  as  well  as  characteristic  slope  or  gradient  of  the  mean 
brightness  distribution.  These  statistics  and  features  of 
one  dimensional  spatial  energy  spectra  are  further  analyzed  to 
determine  rates  of  dissipation  and  diffusion  coefficients  as 
indicated  in  Section  II.  Other  methods  of  calculating  heating 
and  eddy  diffusivity  from  wind  and  temperature  data  are  also 
contained  in  Section  II,  as  is  an  evaluation  of  the  possible 
application  of  video  densitometry  to  extraction  of  the  density 
distributions  from  the  photographs.  Transport  of  the  tracer 
material  by  wind  is  determined  by  triangulation  (which  also 
determines  the  film  plate  scale  or  photogrammetr ic  calibration 
needed  to  determine  turbulence  scales)  as  outlined  in  Section 


III.  Also  described  in  this  Section  is  an  inverse  trian¬ 
gulation  method  utilized  to  recover  trail  altitude  and 
range  information  from  tabulated  wind  profiles  originally 
determined  by  triangulation.  The  task  of  searching  for  a  long 
tern  solar  driven  variation  in  turbopause  altitude  is  discussed 
in  Section  IV. 


SECTION  II 


EDDY  DIFFUSION  AND  HEATING  RATES 

Background 

Turbulent  flow  in  the  atmosphere  is  initiated  and  driven 
by  wind  shears  (Ref  6)  and  by  temperature  instabilities  pro¬ 
duced  by  differential  absorption  of  solar  radiation  by  different 
molecular  and  atomic  species.  Convective  and  turbulent  trans¬ 
port  dominate  over  simple  molecular  diffusion  through  the 
upper  mesosphere  and  into  the  lower  thermosphere,  with  the 
boundary  of  the  turbulent  region  (the  turbopause)  varying 
diurnally  and  seasonally  (and  perhaps  also  over  longer  periods  - 
see  Section  IV)  from  roughly  90  to  120  km.  Computer  codes 
(such  as  Ref  18)  that  model  atmospheric  chemical  species 
concentration  profiles  take  into  account  this  enhanced  mixing 
and  transport  of  species  (that  is,  photochemical  generation 
at  one  altitude  and  transport  to  another  (Ref  19)).  Addition¬ 
ally,  eddy  diffusion  may  affect  the  absolute  atmospheric  density 
and  significantly  affect  vertical  heat  transport. 

Characterization  of  this  turbulent  flow  may  be  accomp¬ 
lished  by  analyzing  surface  brightness  distributions  of 
chemiluminescent  or  sunlight  scattering  tracers  released  from 
rockets.  Fourier  and  configuration  space  statistical  and 
numerical  analysis  methods  are  applied  to  these  surface 
brightness  distributions  to  extract  information  about  the 
turbulent  structure.  A  brief  description  of  the  methods  and 
assumptions  utilized  is  presented  here  (analysis  details  may 
be  found  in  references  1-4). 

The  tracer  material  released  from  the  rocket  (atoms, 
molecules,  or  particulates)  very  quickly  reaches  pressure  and 
temperature  equilibrium  with  ambient  species  and  is  found 
from  experience  to  conform  to  the  flow  field.  It  is  assumed 


Chat  local  fluctuations  in  atmospheric  particle  concentration 
are  traced  by  means  of  this  sunlight  scattering  gas  or  smoke. 

The  spatial  characteristics  of  the  flow  are  frozen  in  time  by 
radiometrically  and  photogrammetrically  calibrated  photography 
from  ground  stations  at  accurately  known  locations.  The  surface 
brightness  distribution  (a  measured  of  line  of  sight  column 
density  of  tracer  serves  as  the  starting  point  for  the  measure¬ 
ments.  Note  that  the  tracer  must  be  optically  thin  to  its  own 
scattered  radiation  for  the  surface  brightness  to  be  a  measure 
of  scale  size  and  tracer  distribution  (usually  assumed  gaussian) 
along  the  line  of  sight. 

Analysis  Method 

Analysis  proceeds  with  microdensitometer  scans  across 
the  images  of  this  frozen  flow  field,  which  are  assumed  (by 
the  ergodic  hypothesis)  statistically  equivalent  to  the  time 
resolved  output  of  a  fixed  physical  sensor  as  the  flow  field 
moves  across  it.  The  resulting  digital  density  information  is 
converted  to  scene  brightness  through  the  H&D  characteristic 
of  the  film,  with  appropriate  corrections  applied  to  account 
for  reciprocity  effects  in  exposure  duration,  spectral  sensi¬ 
tivity  of  the  film,  lens  T-stop,  and  varying  sky  background 
(Ref  4).  Practical  difficulties  encountered  include  film 
grain  noise  limitations  on  spatial  sampling  which  sets  the 
upper  wavenumber  cutoff  for  subsequent  Fourier  analysis;  the 
finite  time  required  to  expose  a  measurable  image  of  the  moving 
cloud  causes  a  smearing  of  the  image  which  again  sets  an  upper 
spectral  limit;  and  the  limited  spatial  extent  of  the  cloud  sets 
an  analogous  cutoff  at  low  wavenumbers  and  also  results  in 
losses  in  statistical  accuracy. 

Further  corrections  are  also  performed  to  account  for  the 
transfer  function  of  the  microdensitometer's  finite  sampling 
aperture  and  for  the  modulation  transfer  function  of  the  micro- 
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dens i tometer  itself,  which  is  a  function  of  the  temporal 
scanning  rate  of  the  instrument.  To  obtain  the  precision 
and  accuracy  needed  for  these  analyses,  the  response  of  the 
data-access ing  instrument  must  be  temporally  and  spatially 
stationary;  for  this  purpose  true  microdensitometer  optics  is 
required,  as  as  explained  in  Reference  11.  (Serious  distortion 
of  the  calculated  power  spectra  can  result  from  nonuniform 
reproduction  of  the  photographically  recorded  data.) 

Photogrammetric  calibration  of  the  images  is  also 
required  (as  described  in  Section  III)  to  determine  the  alti¬ 
tudes  of  the  points  of  measurement  and  to  convert  fiT  alane 
measurements  to  distance  scales  at  the  glow  volume.  me 
sequential  analyses  often  require  a  correction  of  sea  cactor 
due  to  the  range  variation  produced  by  cloud  movement  .  the 
mean  wind  field.  Finally,  the  corrected  brightness  data  are 
analyzed  by  Fourier  (or  other  statistical)  methods.  Fourier 
analysis,  or  calculation  of  one  dimensional  turbulent  energy 
spectra  is  carried  out  by  one  of  three  methods  --  standard 
Fourier  analysis,  "fast"  Fourier  Transform  ( FFT) ,  or  a  simple 
summation  process  (counting  method,  Ref  9)  which  determines 
the  spectral  energy  in  discrete,  nonuniform  wavenumber  bands. 
Each  of  these  methods  requires  large  amounts  of  input  data 
to  achieve  the  desired  statistical  accuracy,  usually  20  to 
50  parallel  microdensitometer  scans  through  the  turbulent 
region  of  interest.  Depending  upon  analysis  method,  the 
average  energy  spectrum  is  obtained  by  segmental  averaging 
or  by  averaging  of  spectra  derived  from  individual  scans. 


Determining  Further  Dynamics 
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exhibits  power  law  dependencies  on  wavenumber,  the  most  well 
known  being  the  Kolmogoroff  or  inertial  spectral  range 
which  implies,  from  dimensional  arguments,  that  turbulent 
energy  is  being  fed  at  a  constant  rate  from  large  cells  into 
small  ones  in  the  subrange  of  cell  dimensions  in  which  this 
slope  obtains.  A  further  application  of  the  energy  spectrum 
is  its  capability  of  isolating  features  of  characteristic 
wavelength  --  for  example  internal  atmospheric  gravity  waves. 

When  an  inertial  spectral  range  is  present  the  heating 
rate  or  rate  of  viscous  dissipation  of  turbulent  kinetic  energy 
c  may  be  found  as  in  Ref  1 .  Further  analysis  may  be  performed 
to  determine  the  kinematic  viscosity  v,  an  eddy  diffusion 
coefficient  K,  and  the  Brunt  Vaisalia  frequency  N2  . 

Considerable  attention  was  given  to  the  persistent 
problem  of  accurate  calculation  of  the  gradient  of  the 
mean  density  structure  an/az,  where  n  is  the  mean  density 
and  z  is  altitude.  The  value  of  an/az  may  be  used  along 
with  the  Fourier  turbuence  spectrum  to  obtain  an  estimate  of 
the  local  heating  rate  from  E^k)  =  a(n/g)  (an/az)  e^/3 ^  -5/3 
where  g  is  gravitational  constant,  k  is  wavenumber,  and  «  is 
a  constant.  Automatic  calculation  of  an/az  from  digital 
microdensitometer  recordings  of  image  film  density  structure 
have  consistently  overestimated  the  gradient  of  the  mean  when 
compared  to  manual  calculation  of  the  gradient  from  two  dimen¬ 
sional  microdensitoneter  recordings.  These  efforts  were 
aimed  at  developing  an  algorithm  which  would  not  overestimate 
the  value  of  the  gradient  of  the  mean  by  typically  factors  of 
20  to  100. 

The  gradient  may  be  manually  determined  by  calculating  the 
slope  (an/az)  of  the  brightness  distribution  (usually 
roughly  gaussian)  transerse  to  the  trail  axis.  Measurements  of 
contour  position  taken  from  the  two  dimensional  microdensitometer 
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plots  of  inage  density  are  plotted  against  corresponding  bright¬ 
ness  data  and  the  slope  calculated  at  the  points  where  the 
surface  brightness  has  fallen  to  1/etfl  of  the  maxinun  bright¬ 
ness.  The  m icrodens itoneter  scanning  aperture  used  to 
produce  these  contour  plots  is  large  (usually  100-200-, n) 
and  the  resulting  plots  are  relatively  immune  to  effects  of 
film  grain  noise  and  small  scale  turbulence  structure  in 
the  image. 

The  automatic  method  uses  digital  microdensitometer 
data  primarily  recorded  for  use  in  determining  Fourier  spectra 
of  turbulence  where  a  small  ('25un)  scanning  aperture  is 
essential  to  allow  calculation  of  the  spectrum  at  high  wave- 
numbers  as  well  as  to  prevent  aliasing  or  "folding"  of  energy 
from  higher  to  lower  wavenumbers.  Thus  small  scale  turbulent 
structure  as  well  as  some  grain  noise  is  included  in  this 
data.  Various  methods  of  filtering  or  smoothing  the  digital 
data  prior  to  calculation  of  the  gradient  were  tested.  Simple 
as  well  as  more  complex  differentiators  were  also  tried  for 
the  actual  calculation  of  the  Dn/z.  Comparisons  to  manual 
calculations  performed  from  equidensity  contour  plots  were 
disappointing,  as  none  of  the  automatic  methods  could  con¬ 
sistently  provide  a  slope  which  compared  to  better  than  a 
factor  of  10.  The  use  of  longer  digital  filters  provided  the 
best  results,  however  filter  length  needed  to  be  several  times 
the  effective  aperture  smoothing  length  used  for  the  equidensity 
plots.  That  is,  filter  length  needed  to  be  on  the  order  of 
one  to  two  millimeters  to  produce  best  results,  which  is  a 
significant  fraction  of  the  total  scan  length  for  these  trans¬ 
verse  scans.  When  parallel  scans  of  digital  data  were  examined 
significant  differences  (factors  of  3  to  5)  in  the  gradient 
were  noted  from  scan  to  scan.  By  averaging  gradients  for  up 
to  50  parallel  digital  scans  automatic  and  manual  reduction 


methods  now  generally  agree  within  acceptable  factors  of  2  to 
3.  Differences  may  now  be  attributed  to  incomplete  smoothing 
of  quantizing  "noise"  in  the  digital  data  and  to  subjective 
smoothing  applied  when  nv  suring  positions  of  contours  on  the 
equidensity  plots. 

Analysis  by  Video  Densitometry 

Although  as  mentioned  earlier,  true  nicrodensitoneter 
optics  is  required  to  obtain  sufficient  accuracy  for  Fourier 
analysis  for  turbulence,  the  slow  data  acquisition  rate  of  the 
system  used  led  to  an  investigation  of  the  possible  use  of  an 
existing  computerized  video  based  densitometer.  The  spatial 
sampling  characteristics  (i.e.  sampling  increment  uniformity) 
of  the  system  were  found  to  be  adequate  for  this  task.  However 
the  spatial  variation  of  density  response  for  this  full  field 
illumination  instrument  showed  both  correctable  and  apparently 
uncorrectable  aberrations  in  its  measurement  capability. 

Vignetting  (shading)  due  to  the  densitometer  camera 
lens,  nonuniformity  of  the  illumination  source,  and  also  the 
nonuniformity  of  response  over  the  video  photocathode  are  all 
correctable.  In  fact,  a  single  correction  step  (pixel  by 
pixel)  using  a  previously  recorded  reference  field  is  all  that 
is  needed  to  effect  this  time-stationary  correction.  The 
reference  field  is  simply  the  image  generated  by  the  video 
system  with  no  film  in  its  field  of  view. 

The  virtually-uncorrec table  problem,  which  is  the  usual 
difficulty  associated  with  full  field  illumination  --  the 
contribution  of  scattered  light  --  is  that  the  density 
response  of  the  system  at  each  pixel  is  highly  dependent  upon 
the  micro-density  (or  light  absorption)  of  the  surrounding 
pixels.  It  is  conceivable  that  for  structurally  very  simple 
backgrounds  and  objects  a  correction  for  this  flare  light 


background  could  be  determined.  However  images  of  turbulence 
are  certainly  not  in  this  category  and  it  is  unlikely  that 
suitable  algorithms  for  correction  could  be  developed.  New 
camera  sensor  techniques  are  not  expected  to  eliminate  this 
problem.  (It  remained  when  a  CCD  array  camera  was  substituted 
for  the  vidicon).  Even  a  linear  detector  array,  which  mechan¬ 
ically  scans  in  one  direction,  would  be  susceptible  to  the 
flare  light  effect.  Flare  is  controlled  in  a  true  microdensi- 
tometer  by  use  of  narrow  field  illumination,  and  a  preslit  to 
further  limit  the  illuminated  image  area  to  only  that  micro¬ 
area  viewed  by  the  narrow  field  objective  lens. 

A  further  nonuniformity,  which  involved  temporal  non- 
stationarity  of  density  response  to  individual  images,  was 
noted.  The  measured  optical  density  would  vary  cyclically, 
apparently  in  response  to  changes  in  ambient  room  temperature. 
Although  relatively  severe  (30-40%  variation  at  times)  ,  this 
could  be  corrected  by  stabilizing  the  room  temperature. 

While  most  of  the  problems  associated  with  the  use  of 
video  based  microdensitometry  can  be  adequately  corrected, 
flare  effects  due  to  the  combination  of  full  field  illumination 
scattering  within  the  sample,  and  acceptance  angle  of  the 
sensor  are  apparently  uncorrectable .  Although  such  systems 
provide  extremely  rapid  digitization  of  photographically 
recorded  data  the  spatial  non  stationarity  of  photometric 
response  is  considered  unacceptable. 


SECTION  III 


PHOTOGRAPHIC  TRIANGULATION 

|  Determinat ion  of  turbulence  parameters  from  photographic 

records  of  atmospheric  tracer  releases  requires  an  accurate 
knowledge  of  the  true  scale  sizes  at  the  glow.  This  is  usually 
accomplished  by  performing  photogrammetr ic  triangulation  from 
|  two  or  more  photographic  projections  to  determine  its  three- 

dimensional  spatial  location,  and  hence  its  range  from  the 
camera  which  defines  the  conversion  from  image  scale  to 
dimensions  at  the  object.  The  methods  used  to  perform  these 
I  scale  calibrations  are  detailed  in  References  12  and  13,  with 

further  information  on  their  applications  in  References  14  and 
15.  A  brief  description  of  the  technique  is  presented  here. 

Photographic  triangulation  to  distant  objects  is  based 
upon  accurate  knowledge  of  the  location  and  orientation  of 
each  camera  photographing  the  object.  The  location  is  usually 
expressed  in  latitude,  longitude,  and  altitude  (referenced  to 
an  earth  ellispoid  model) ,  which  are  determined  by  surveying 
the  camera  positions.  Positions  may  sometimes  be  scaled  from 
USGS  topographic  maps,  although  the  lower  accuracy  may  severely 
degrade  the  triangulation  results.  The  orientation  of  the 
camera  (the  angles  which  describe  the  viewing  direction)  is 
determined  by  analysis  of  images  which  have  recorded  the  star 
field  at  accurately  known  times.  The  film  frames  with  star 
images  are  semiautomatically  aligned  on  the  film  digitizer 
(Ref  14)  such  that  an  orthogonal  coordinate  system  defined 
by  the  film  fiducial  markers  (coordinate  (0,0)  is  the  camera 
optic  axis)  coincides  with  the  digitizer  axes.  A  digital  star 
catalog  is  searched  to  determine  which  stars  are  in  the  field 
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of  view  and  the  digitizer  semiautomatically  locates  each  star 
image  and  records  its  film  plane  coordinates  along  with  equa¬ 
torial  coordinates  (right  ascension  and  declination)  of  the 
star  corrected  to  the  precise  time  of  the  exposure.  The  film 
plane  and  equatorial  coordinates  for  up  to  20  stars  are 
related  through  vector  equations  which  are  solved  to  deter¬ 
mine  the  lens  focal  length  (using  star  triplets)  and  also  the 
direction  in  space  of  the  camera  optic  axis  (using  pairs  of 
stars).  The  orientation  of  the  optic  axis  is  calculated  in 
right  ascension  and  declination,  and  also  an  angular  rotation 
needed  to  bring  the  star  and  image  locations  into  coincidence. 
It  is  often  more  convenient  either  for  ease  of  visualization 
or  for  further  calculation,  to  express  (the  orientation  in 
some  other  coordinate  system  (such  a  horizon  or  geocentric) ; 
this  is  easily  accomplished  through  coordinate  transformation 
of  the  optic  axis  line  of  sight  vector  (rotations  about  co¬ 
ordinate  axes)  . 

Time-coincident  frames  from  two  photographic  observation 
sites  are  digitized  using  the  same  automatic  scanning  densi¬ 
tometer  as  used  for  star  frames.  This  system  locates  and 
digitizes  the  film  plane  coordinates  of  the  centerline  of  the 
glow  image  for  use  by  the  triangulation  routines.  Turbulent 
regions  of  releases  often  do  not  have  readily  identifiable 
centerlines  and  automatic  reduction  of  the  film  density 
distribution  transverse  to  the  trail  axis  to  locate  the 
center  is  ineffective.  They  may  however  have  "blobs"  of 
tracer  material  which  are  uniquely  identifiable  from  each 
site  whose  coordinates  may  be  digitized.  Further,  when 
only  approximate  altitude  information  is  needed  (to  associate 
altitudes  with  energy  spectra  or  other  turbulence  uaraneters) 
a  series  of  small  "dots"  or  marks  may  be  subjectively  placed 
on  each  film  (or  an  overlay)  along  the  trail  "centerline." 
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While  positions  of  these  dots  are  easily  digitized,  the 
resulting  triangulation  is  imprecise  (dots  on  one  image  will 
generally  not  correspond  closely  to  dots  on  the  other  triangu¬ 
lation  image)  and  interpolation  is  required  to  associate  space 
positions  with  each  dot  location. 

The  digital  film  plane  positions  define  a  series  of 
vectors  to  the  tracer  glow  from  each  photographic  site  and 
are  used  to  define  the  spatial  location  of  the  glow.  An 
iterative  procedure  is  used  to  determine  the  best  match  of 
each  site  1  line-of-s ight  vector  with  a  line-of-sight  vector 
from  site  2.  The  point  of  intersection  (or  closest  approach) 
of  each  vector  pair  is  found  in  the  geocentric  coordinate 
system,  converted  to  latitude,  longitude,  altitude  in  the 
geographic  system,  and  tabulated  for  association  with  a  par¬ 
ticular  point  on  each  film  image.  Although  conceptually 
the  technique  for  triangulation  is  straightforward  there  are 
problems  associated  with  the  implementation  of  the  method, 
some  of  which  are  detailed  in  Ref  15.  These  problems  are 
primarily  with  the  design  of  the  iterative  routines  which 
determine  the  "best  match"  pairing  of  points  using  a  dihedral 
angular  matching  criterion  and  with  resolving  trail  segments 
which  lie  nearly  along  a  uniform  line  of  sight  when  viewed 
from  one  of  the  photographic  projections. 

Triangulation  is  also  the  first  step  in  remote  measure¬ 
ment  or  the  wind  field  at  high  altitudes.  The  trail  position 
is  first  determined  for  several  photographic  frame  pairs 
extending  over  total  durations  of  one  to  two  minutes.  Each 
altitude  profile  of  position  is  spline  interpolated  to  a 
convenient  equal  altitude  grid  and  least  square  analysis  of 
trail  position  as  a  function  of  time  is  applied  to  each 
altitude  to  yield  average  positional  changes,  average  velocity 
components,  and  their  associated  errors. 


An  inversion  of  the  triangulation  method  has  been  used 
as  a  reconstruction  technique  to  allow  identification  of 
altitudes  on  films  of  previously-made  releases  for  which 
much  of  the  supporting  information  has  been  lost  or  discarded, 
or  for  which  only  one  of  the  original  triangulation  films  is 
available.  These  chemical  releases  had  been  used  primarily  to 
determine  molecular  diffusion  coefficients  and  wind  field 
(using  overlays  with  manually  placed  dots).  When  the  overlays 
or  the  computer  printout  of  dot  number  and  its  triangulated 
altitude  are  missing  it  is  still  possible  to  use  the  profile 
of  wind  components  (which  are  available  in  the  analyzed  data 
reported)  to  assign  altitudes  along  the  trail  image. 

Utilizing  the  known  trajectory  of  the  release  rocket 
and  the  known  camera  parameters  (focal  length,  azimuth,  eleva¬ 
tion,  altitude,  latitude  and  longitude)  a  computer  program 
first  constructs  a  model  release  along  the  trajectory.  This 
line  release  is  then  transported  using  the  known  wind  profile 
and  at  predetermined  increments  of  time  (matched  to  the  photo¬ 
graphic  times)  vector  triangulation  determines  the  projection 
of  the  release  into  the  film  plane  of  the  particular  camera. 
Computer  generated  plots  of  the  image  plane  at  1  km  altitude 
intervals  were  then  used  as  overlays  to  establish  an  altitude 
grid  on  the  release  film  (Fig  1).  The  major  sources  of  error 
in  this  inverse  triangulation  arise  from  uncertainty  in  site 
locations  and  the  generally  few  points  on  the  trajectory  which 
are  known  (often  these  points  used  are  not  actually  on  the 
trajectory).  Estimated  altitude  errors  have  been  in  the  1-4  km 
range  when  using  this  method. 
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SECTION  IV 


TURBOPAUSE  HEIGHT  ANALYSIS 

Analyses  were  performed  in  an  effort  to  characterize  the 
height  of  the  turbopause  as  a  function  of  several  geophysical 
variables.  The  principal  purpose  of  this  analysis  was  to 
investigate  long  term  (-16  years)  variation  of  turbopause 
height  ( Zt)  as  a  function  of  the  10.7  cm  radiation  from  the 
sun  (F10.7),  and  the  planetary  magnetic  disturbance  indices 
Kp  and  Ap.  It  has  been  suggested  (Ref  16)  that  there  is  a 
strong  solar  cycle  variability  of  Zj  associated  with  long  term 
changes  in  the  10.7  cm  flux  and  possibly  also  associated  with 
the  Kp  and  Ap  indices,  which  themselves  also  correlate  with 
the  11  year  solar  cycle. 

The  data  base  used  for  analysis  is  comprised  of  measure¬ 
ments  of  turbopause  height  determined  from  AFGL  chemical 
release  programs  conducted  at  Eglin  Air  Force  Base  from  1962 
to  1976  and  from  heights  deducted  from  Ar/N2  ratios  measured 
by  rocketborne  mass  spectrometers  at  Heiss  Island  from  1968  to 
1977.  Associated  with  each  Zj  is  the  Fio.7  flux,  Kp  and  Ap 
indices . 

The  most  desirable  method  of  analysis  for  extraction 
of  dependencies  of  Zj  on  many  variables  is  multivariate  re¬ 
gression  analysis.  Such  an  analysis  would  determine  how 
well  (in  terms  of  regression  coefficients  and  standard 
deviations)  the  measured  data  corresponds  to  hypothesized 
functional  forms  for  the  variation  of  each  independent 
parameter.  This  was  not  employed  because  of  the  limited 
number  of  measurements  relative  to  the  number  of  variables. 
Confidence  in  the  results  would  be  greatly  diminished 
because  of  large  standard  deviations. 


Our  anproach  to  determining  any  long  term  variations  was  to 
attempt  to  remove  the  sources  of  short  term  variation  (semi 
diurnal,  diurnal,  seasonal)  from  the  data  and  then  examine  the 
remaining  trends,  if  any.  The  data  were  hence  divided  into 
winter  and  summer  groupings  (16  Sept  to  15  Mar  and  16  Mar  to 
15  Sept  respectively)  similar  to  References  16  and  17,  and  an 
approximate  form  for  the  diurnal  variation  used  to  modify  the 
observed  Z  values  such  that  all  would  be  representative  of  the 
height  at  a  uniform  local  time  (0600  for  winter,  0300  for 
summer).  Each  of  these  times  correspond  to  an  altitude  of 
105  km.  The  forms  for  diurnal  variation  were  semiquantitat ively 
chosen  since  the  majority  of  the  data  is  clustered  near  0300 
hours  and  2000  hours.  The  preponderance  of  measurements  near 
twilight  results  from  conditions  necessary  to  photograph  high 
altitude  sunlight  scattering  releases;  good  contrast  may  only 
be  obtained  when  the  observation  camera  is  in  darkness  while 
the  tracer  material  is  still  sunlit.  Some  of  the  releases 
were  chemi  or  self  luminescent  allowing  nighttime  measurements 
to  be  performed.  Semi  diurnal  corrections  were  not  attempted 
due  to  this  clustered  nature  of  the  measurements.  Furthermore, 
after  diurnal  corrections  were  applied  to  the  seasonally 
separated  Z-p  values,  multiple  measurements  performed  within  the 
same  year  had  their  Z-p  and  Fio.7  values  averaged. 

Figures  2  and  3  show  the  variation  of  Z-p  as  a  function 
of  local  time  for  the  summer  and  winter  groups  respectively. 

Each  figure  also  shows  the  diurnal  form  which  was  used  for 
correction  purposes.  It  is,  of  course  possible  that  some 
of  the  apparent  diurnal  change  could  due  to  long  term  effects 
associated  with  the  limited  number  of  samples.  The  nighttime 
rise  of  the  turbopause  is,  however,  clearly  demonstrated  by 
the  series  of  chemical  release  measurements  (dashed  line- 
winter)  made  on  a  single  night  and  hence  only  the  specific 
form  of  the  diurnal  variation  may  be  questioned. 


Figure  2.  Turbopause  height  data  for  the  summer  grouping  with  diurnal  variation 
curve  estimated.  Dots  represent  data  measured  at  Eglin  AFB  and 
squares  at  Heiss  Island. 


Figure  3.  Turbopause  height  data  for  the  winter  grouping  with  diurnal  variation  curve 

estimated.  Dashed  curve  shows  the  nighttime  rise  determined  from  experiments 
-  —  nacXomad  night  C1962).  Dots  represent  data  measured  at  Eglin  AFB 


Figure  4.  Turbopause  height  as  a  function  of  F10<?,  summer.  Data  have  been 
corrected  for  diurnal  variation  and  averaged  by  year. 


The  behavior  of  seasonal  turbopause  height  by  year 
as  a  function  of  Fio.7  is  shown  in  Figures  4  and  5.  These 
data  have  been  corrected  for  diurnal  variation  and  Z-p  and 
flux  values  averaged  by  year.  These  data  with  year  associated 
with  each  measurement  appear  to  show  a  hysteresis  or  lagging 
of  the  altitude  response  to  changes  in  the  decimeter  flux. 

This  is  most  clearly  demonstrated  by  the  clockwise  circulation 
in  the  summer  data;  winter  data  appear  to  circulate  counter¬ 
clockwise  . 

Results  of  further  attempts  to  remove  uncertainty  due  to 
latitude  and  possible  short  term  changes  in  flux  are  shown  in 
Figures  6  and  7.  An  obvious  systematic  altitude  differential 
between  the  30°N  and  80°N  data  was  compensated  by  estimating 
the  average  altitude  shift  needed  to  make  the  diurnal  variation 
curves  representative  of  either  the  30°N  or  80°N  data.  Further¬ 
more  it  was  assumed  that  altitude  changes  in  response  to  changes 
in  F  i  o .  7  d°  not  occur  rapidly  and  a  smoothed  average  flux 
was  substituted  for  the  yearly  seasonal  averages.  Note  that 
the  direction  of  circulation  is  now  approximately  clockwise 
for  both  summer  and  winter  data.  The  standard  deviations  of 
the  mean  altitude  are  2.5  km  and  7  km  for  winter  and  summer 
respectively.  Both  sets  of  data  have  excursions  barely  exceed¬ 
ing  ±  1  sigma,  and  while  they  may  be  suggestive  of  a  long 
term  dependence  upon  F-|o.7*  they  also  may  be  considered  as 
samples  of  a  normal  distribution  of  altitudes  (t  test  at  95% 
confidence).  Further  analysis  could  be  conducted  vising  this 
data,  especially  the  inclusion  of  the  Kp  and  Ap  values  which 
presumably  could  reinforce  the  result  by  demonstrating  similar 
Z'p  dependence  upon  other  solar  derived  parameters. 


Summer  turbopause  height  corrected  for  observation  latitude  as 
a  function  of  smoothed  10.7  cm  solar  flux. 


SECTION  V 


CONCLUSIONS 

Techniques  developed  and  refined  over  the  past  several 
years  have  been  applied  to  extract  and  analyze  data  from 
photographic  records  of  chemical  releases  to  provide  turbulent 
diffusion  and  dissipation  in  the  region  of  the  turbopause. 
Fourier  turbulence  spectra,  the  gradient  of  mean  density 
distribution  an/dz,  the  rate  of  change  of  turbulent 
density  fluctuations,  and  horizontal  wind  and  shears  were 
utilized  to  determine  diffusion  and  heating. 

Since  the  gradient  of  the  mean  brightness  distribution 
when  combined  with  numerical  parameters  of  the  energy  spectra 
has  provided  the  most  reliable  measurements  of  these  quanti¬ 
ties,  methods  to  refine  the  calculation  of  gradient  were 
investigated.  Several  techniques  were  applied  to  the  tasks 
of  automatically  computing  the  gradient  from  the  digital 
photographic  records,  none  of  which  would  consistently  yield 
results  strictly  comparable  to  manually  derived  values. 
Remaining  factors  of  2  to  3  difference  are  probably  attribut¬ 
able  to  insufficient  smoothing  of  quantizing  "noise"  of  the 
digital  data  and  to  subjective  smoothing  done  during  manual 
reduction  of  the  gradients.  Two  dimensional  spatial  filtering 
of  the  digital  data  may  provide  for  removal  of  this  "noise," 
or  an  interactive  approach  could  be  utilized  to  simulate  the 
manual  reduction  method. 

The  search  for  a  more  rapid  way  to  extract  density 
information  from  the  photographs  led  to  examination  of  video 
access  methods.  Although  video  techniques  are  suitable  for 
some  applications  which  require  quantitative  measurement  of 
film  density,  currently  available  corrections  for  "crosstalk" 


are  not  adequate  to  provide  suitably  accurate  turbulence 
measurements.  Primarily,  flare  light  effects  due  to  full 
field  illumination  vary  with  the  densities  and  spatial  content 
of  the  image,  making  calibration  difficult,  if  not  impossible, 
as  image  complexity  increases. 

Identification  of  the  altitude  (as  well  as  establishing 
the  film  plate  scale)  where  turbulence  measurements  were  made 
was  done  by  triangulation  when  sufficiently  accurate  records 
of  site  locations  and  other  release  data  were  available. 
Alternatively,  a  less  accurate  reconstruction  technique 
(inverse  triangulation)  was  employed  when  adequate  calibration 
data  were  not  available. 

Analysis  was  performed  to  establish  a  probable  long 
term  relationship  between  turbopause  height  and  solar  10.7 
centimeter  flux.  The  results  suggest  a  lagging  11  year  cycle 
of  motion  of  the  turbopause  altitude  in  response  to  the  Fio.7 
flux.  Confidence  in  the  results  is  not  great  because  of  the 
limited  data  available.  Analysis  of  the  corresponding  Kp  and 
Ap  magnetic  index  data  could  possibly  establish  a  similar 
solar  activity  altitude  dependence. 

The  results  of  the  data  analysis  for  turbulence  para¬ 
meters  were  reported  to  AFGL  scientists  as  they  were  obtained, 
and  are  on  file  at  the  LYT  Branch  of  AFGL. 
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